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Abstract 

We consider resonant inverse Compton scattering of thermal photons 
by secondary particles above the pulsar polar gap. At neutron star tem- 
peratures higher than 10 5 K the process appears to be an essential en- 
ergy loss mechanism for the particles. The distribution function of the 
secondary plasma particles is found to be strongly affected by the scat- 
tering. It becomes two-humped implying the development of two-stream 
instability. The resonantly upscattered Compton photons are found to 
gain energy of 1-10 MeV forming an additional component in the pulsar 
gamma-ray spectrum. The corresponding gamma-ray flux is estimated as 
well. 

1 Introduction 

Rotation of a highly magnetized neutron star is known to induce a strong electric 
field, which intensely accelerates charged particles. According to the customary 
polar gap models (Ruderman & Sutherland 1975, Arons & Scharlemann 1979), 
the acceleration takes place near the neutron star surface above the polar cap, 
the Lorentz-factor of the primary particles increasing up to ~ 10 6 . The parti- 
cles move along the magnetic lines and emit curvature photons, which initiate 
pair-production cascade. The first electron-positron pairs created screen the ac- 
celerating electric field, so that at higher altitudes the particle energy remains 
unaltered; the typical Lorentz- factors of the secondary plasma are ~ 10 — 10 4 . 

Recent observations testify to thermal soft X-ray emission from some pulsars 
indicating that the neutron stars can be rather hot, T ~ 5 • 10 5 K, while the 
polar caps can have still higher temperatures scaling a few times 10 6 K (Cordova 
et al. 1989, Ogelman 1991, Halpern & Holt 1992, Finley et al. 1992, Halpern 
& Ruderman 1993, Ogelman & Finley 1993, Ogelman et al. 1993, Yancopoulos 
et al. 1994, Ogelman 1995, Greiveldinger et al. 1996). Such high temperatures 
of the neutron star surface are also predicted theoretically (Alpar et al. 1984, 
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Shibazaki & Lamb 1989, Van Riper 1991, Page & Applegate 1992, Umeda et 
al. 1993, Halpern & Rudcrman 1993). The thermal X-ray photons should suf- 
fer inverse Compton scattering off the primary particles in the polar gap. In 
the presence of a strong magnetic field the scattering cross-section is essentially 
enhanced, if the photon energy in the particle rest frame equals the cyclotron 
energy (Herold 1979, Xia et al. 1985). For the pulsars with hot polar caps the 
resonant Compton scattering in the polar gap was found to be rather efficient 
(Kardashev et al. 1984, Xia et al. 1985, Daugherty & Harding 1989, Dermer 
1990, Sturner 1995, Chang 1995). Firstly, it was recognized as an essential 
mechanism for energy loss of primary particles accelerating in the polar gap. 
Secondly, inverse Compton scattering was found to condition the gap forma- 
tion, since the pair production avalanche may be triggered by the upscattered 
Compton photons rather than by curvature photons (Zhang & Qiao 1996, Qiao 
& Zhang 1996, Luo 1996, Zhang et al. 1997). 

As shown by Sturner (1995), given typical values of neutron star temperature 
and surface magnetic field, the resonant Compton scattering is the strongest for 
particle Lorentz-factors ~ 10 2 — 10 3 . Theoretical models (Van Riper 1991, Page 
& Applegate 1992, Umeda et al. 1993) suggest that the typical temperatures 
of the neutron star surface are as high as a few times 10 5 K. Hence, the scat- 
tering is likely to be essential for the secondary plasma in most of the pulsars. 
Daugherty & Harding (1989), Zhang et al. (1997) traced the evolution of the 
Lorentz-factor of secondary particles on account of magnetic inverse Compton 
scattering above the polar gap. Our aim is to investigate in more detail the in- 
fluence of resonant inverse Compton scattering on the parameters of secondary 
plasma. For the resonant character of the scattering, the evolution of particle 
Lorentz-factor with the distance depends strongly on the initial particle energy. 
Since the distribution function of the secondary plasma is generally believed 
to be rather broad (the Lorentz-factor ranges from 10 to 10 4 ), its evolution 
on account of the Compton scattering is of a great interest. It will be shown 
that only particles with Lorentz-factors between 100 and 1000 are essentially 
decelerated in the course of the resonant scattering forming a sharp peak at low 
energies. Particles with larger Lorentz-factors are not decelerated at all. Thus, 
the resultant distribution function of secondary particles becomes two-humped, 
giving rise to the two-stream instability. 

In Sect. 2 we examine how the resonant inverse Compton scattering affects 
the Lorentz-factor of a secondary particle at various pulsar parameters. Sec- 
tion 3 is devoted to studying the evolution of the distribution function. The 
conditions for the development of the two-stream instability are also discussed. 
In Sect. 4 we estimate the gamma-ray luminosity caused by the upscattered 
Compton photons. Section 5 contains a brief summary. 
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2 Deceleration of secondary particles due to res- 
onant inverse Compton scattering 

Consider the flow of secondary plasma streaming along the open magnetic field 
lines above the polar gap. The neutron star is supposed to emit a blackbody 
radiation, which is scattered by the plasma particles. In a strong magnetic field 
the scattering is particularly efficient, if the photon energy, emc 2 , satisfies the 
resonance condition: 

£ 7 (l-/3cos0) =e B , (1) 

where 7 is the Lorentz-factor of the scattering particles, (3 the particle velocity 
in units of c, the angle the photon makes with the particle velocity, e = B/B cr , 
with B cr = m 2 c 3 /he — 4.414 • I0 13 G. At the distance z from the neutron star 
the rate of energy loss due to the resonant inverse Compton scattering can be 
written as (Sturner 1995): 
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Here Tq is the neutron star temperature in units of 10 6 K, Byi{z) the magnetic 
field strength in units of 10 12 G, 9 c (z) the maximum incident angle of photons, 



where R is the neutron star radius. Provided that the magnetic field is dipolar, 
B 12 (z) oc (l + z/R)- 3 . 

Numerical solutions of Eq. (2) are presented in Fig. 1. 

One can see that the resonant inverse Compton scattering is significant up 
to z ~ R and it affects the particle Lorentz-factor essentially. Note that the 
curves for different initial Lorentz-factors are not similar to each other. In 
agreement with Eq. (1), the larger the particle Lorentz-factor the lower is 
the energy of resonantly scattered photons. At 7 o = 3000 the energy of the 
resonant photons is below e max — 2.82fcT/(mc 2 ), which corresponds to the 
maximum in the photon distribution. As 7 is decreasing with the altitude, the 
resonant energy increases. Hence, the photon spectral density increases and the 
scattering becomes more efficient. However, at distances z ~ R it ceases due to 
the essential shrinking of the solid angle subtended by the neutron star. If the 
initial Lorentz-factor is 100, the energy of resonant photons is above e max and 
further increases with the altitude, so that the scattering gradually ceases. 

Theoretical models predict that the polar cap region should be significantly 
hotter than the rest surface of the neutron star {e.g. Cheng & Ruderman 1980, 
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Arons 1981, Alpar et al. 1984, Umcda et al. 1993, Luo 1996). However, the 
luminosities observed from the hot spots appear to be too small (Finley et 
al. 1992, Halpern & Ruderman 1993, Becker & Trumpcr 1993, Yancopoulos 
et al. 1994, Greiveldinger et al. 1996). The latter is usually interpreted as a 
consequence of small hot spot radii (e.g. Yancopoulos et al. 1994, Greiveldinger 
et al. 1996). As is evident from Fig.l, for typical hot spot parameters the 
particle energy loss is mainly determined by the scattering of the photons from 
the whole neutron star rather than by the scattering of hot spot photons. So 
hereafter we take into account only the photons from the whole neutron star 
surface keeping in mind that the influence of hot spot photons can somehow 
alter our quantitative results, while the qualitative picture should remain the 
same. 

It should be pointed out that the evolution of Lorentz-factor with the al- 
titude shown in Fig. 1 is somewhat different from that reported by Zhang et 
o/.(1997)(see their Fig. 1). These authors claim that at distances z ~ fewi? the 
particles suffer severe nonresonant magnetized scattering, which leads to the 
drastic decrease of the final Lorentz-factor. In fact, the rate of energy loss on 
account of nonresonant magnetic scattering increases with decreasing the field 
strength as B^ 2 (z) (see Eq. (5) in Sturner 1995). However, this equation is 
applicable only if in the particle rest frame the cyclotron energy exceeds the 
energy of most of photons, Eb > £7(1 — (3cos8 c ). For the photon energies at 
the peak of the Planck distribution, e w 2.82fcT, the latter condition can be 
rewritten as 4 • 1CT 2 B 12 (z) / '(T 6 7 3 (l - /3cos6» c )) > 1, with 73 = 7/IO 3 . Taking 
Tg = 0.5, 73 = 3, B12 = 0.1 (Zhang et al. 1997, Fig. la), one can obtain that 
at z = lOi? the left-hand side of this inequality is 5 • 10~ 5 , while at z = it 
equals 3 • 10~ 4 . So even at the stellar surface most of the photons scatter in the 
Thomson regime rather than in the magnetic one. The rate of energy loss out 
of the Thomson scattering is found to be 

^ = -3O7T 6 4 (l-/?cos0 c ) 3 s-\ (4) 

showing that this process is inefficient. Thus the resonant inverse Compton scat- 
tering is the only significant energy loss mechanism for the secondary particles 
in pulsar magnetosphcres. 

We believe that the particles start deceleration just above the polar gap. In 
general the gap thickness is supposed to be of the order of the polar cap radius 
(Ruderman & Sutherland 1975) or even larger (Arons & Scharlemann 1979). 
The energy loss of secondary particles due to the resonant inverse Compton 
scattering above the gap should certainly be influenced by the adopted value 
of the gap height. In Fig. 2 we show the final Lorentz-factors versus the gap 
height. One can see that the dependence is sufficiently weak, therefore, below 
we fix z at 10~ 2 i?. 

In contrast with the gap height, such pulsar parameters as the star temper- 
ature and magnetic field strength can influence particle deceleration essentially. 



4 



In Fig. 3 we present the dependences of the normalized final Lorentz-factor on 
the neutron star temperature at various magnetic field strengths. Apparently, 
at the temperatures < I0 5 K the scattering is inefficient yet, while at higher 
temperatures it becomes significant. Note that at various initial Lorentz-factors 
(70 = 100 and 3000) the scattering efficiency versus magnetic field strength is 
essentially different. The particles with 70 = 100 resonantly scatter the photons 
from the Wien region. Then the weaker the field strength, the lower is the en- 
ergy of resonant photons and, correspondingly, the higher is the photon spectral 
density and the larger is the particle energy loss (see Fig. 3a). At 70 = 3000 
the resonant energy lies in the Raylcigh- Jeans region. So the scattering is more 
efficient for higher magnetic field strengths, since the photon spectral density 
increases with the energy (see Fig. 3b). 

For the resonant character of the scattering the particle energy loss depends 
strongly on the initial Lorentz-factor. Figure 4 shows the final Lorentz-factor 
versus the initial one for various neutron star temperatures and magnetic field 
strengths. At 70 ~ 10 as well as at 70 ~ 10 4 the resonant scattering appears to 
be inefficient. Provided that the scattering is intense (70 ~ I0 2 — I0 3 ), the final 
Lorentz-factor appears to be independent of the initial one. According to Fig. 
4a, the length of the plateau increases with the star temperature. In fact, the 
higher temperature implies the larger amount of photons at every energy, the 
scattering becoming more efficient. As can be seen from Fig. 4b, the increase 
of the magnetic field strength leads to the shift of the plateau toward a higher 
energy; this is certainly consistent with the resonance condition (1). 



3 The distribution function of the secondary plasma 
as a result of resonant inverse Compton scat- 



Since the energy loss depends essentially on the initial particle energy, we are to 
investigate the evolution of the distribution function of secondary particles as a 
result of the resonant inverse Compton scattering. Conservation of the number 
of particles along a phase trajectory implies that 



where f(z, 7) is the particle distribution function and the subscript "0" refers 
to the initial values. So looking for the phase trajectories of individual particles 
one can reconstruct the distribution function at any height z. We are particu- 
larly interested in the final distribution function arising at distances, where the 
resonant scattering ceases. 

Let us begin with the evolution of a waterbag distribution function: 
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with 7 m i n = 10, 7™^ = 10 . The final distribution functions at various pulsar 
parameters are plotted in Fig. 5. The particles with 7 ~ 10 2 — 10 3 are essentially 
decelerated due to the scattering, the final energies becoming equal (see also Fig. 
4). These particles form the sharp peak at 7 < 10 2 . For the Lorentz-factors of 
a few thousand the scattering becomes inefficient and the distribution function 
remains almost unaltered. Thus, the resonant inverse Compton scattering leads 
to the two-humped distribution function of the secondary plasma. At higher 
neutron star temperatures the main peak of the function shifts towards the lower 
energies and the humps become more prominent. The magnetic field strength 
variation also results in the shift of the main peak. For the more realistic 
distribution function resembling that found by Arons (1980), 

exp(-lO^), 10< 7 <7m, 
f(l)={ (7/7m)- 3/2 , 7™<7<7c (5) 

(lch m r 3/2 exp (-^) , 7c < 7 < 10 4 , 

with 7 m = 10 2 , 7 C = 10 3 - 5 , the evolution on account of the resonant inverse 
Compton scattering is qualitatively the same (see Fig. 6). 

We next examine the dispersion properties of the plasma with the evolved 
distribution function. For simplicity let us assume that the plasma consists of 
the two particle flows characterized by the number densities n a , rib and by the 
Lorentz-factors -y a , 7^ (j a < 7;,). Given the infinitely strong magnetic field, the 
dispersion relation is as follows (see, e.g. Lyubarskii 1995): 



, ,2 , ,2 



= I- (6) 



Here ui is the frequency, k the wave number, v a ,b are the particle velocities, Lo pa ,b 
the plasma frequencies given by the customary expression: 



m 



where e is the electron charge, m the electron mass. As is evident from Eq. (6), 
the dispersion properties of the plasma are mainly determined by one of the 
particle flows on condition that 

-z » ^ W 

7a 7 6 

rather than n a 3> rib. This is because of the great inertia of the fast particles 
performing one-dimensional motion in the superstrong magnetic field. For the 
distribution functions plotted in Fig. 5 and 6 n a and n arc comparable, while 
7a/76 ~ 10 _1 — 10~ 2 . Therefore the low-energy particles of the main peak 
almost completely determine the dispersion properties of the plasma. 
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The two-humped distribution function implies the possibility of the two- 
stream instability. Provided that the contribution of one of the plasma flows 
to the plasma dispersion is small (i.e. Eq. (8) is valid), the growth rate of the 
instability takes the form (Lominadze & Mikhailovskii 1979, Cheng & Ruderman 
1980): 

V3 ,, 



m*,« ^1 (9) 

The two-stream instability results in an essential development of initial pertur- 
bations on condition that 

-Imw > 10, (10) 

c 

where R c is the characteristic scale length for the increase of perturbations. 

It is convenient to normalize the number density of the secondary plasma by 
the Goldreich-Julian charge density: 

where k is the multiplicity factor of the secondary plasma, P the pulsar period. 
The latter equation can be rewritten as: 

n = 6.25 • 10 13 P- 1 K 3 -Bi2(l + z/R)~ 3 cnT 3 , (12) 

where K3 = k/10 3 . Using Eqs. (9) and (12) in Eq. (10) we reduce the condition 
for the efficient instability development to the form: 

H:dh} >4.10- 4 . (13) 



PRc7% 3 la2 



Here R C 7 = i? c /10 7 cm, 7^3 = 7;,/10 3 , -f a 2 = 7a/10 2 and it is assumed that 
(nb/ria) 1 / 3 « 1. As can be seen from Eq. (13), at typical pulsar parameters 
the two-stream instability can develop readily providing the increase of plasma 
oscillations. The latter can be transformed into electromagnetic waves, thus 
giving rise to pulsar radio emission. 

It should be mentioned that two-steram instability has always been one of 
the most popular mechanisms for pulsar radio emission. For more than two 
decades a number of scenarios for the instability development were proposed. 
The first and the most natural one involves the instability caused by the Hows 
of primary and secondary pulsar plasma (Ruderman & Sutherland 1975). How- 
ever, the growth rate of this instability appears to be too small because of 
enormous inertia of the high-energy primary particles (Benford & Buschauer 
1977). Cheng & Ruderman (1977) considered the two-stream instability aris- 
ing in the secondary plasma due to the difference in the velocities of electrons 
and positrons moving along the curved magnetic lines. However, this difference 
is insufficient to cause the instability, since the particle distribution functions 



7 



are too broad (Buschauer & Benford 1977). Lyubarskii (1993) suggested that 
current and charge density adjustment in pulsar magnetosphere leads to the 
backward particle flow, which causes intense two-stream instability. However, 
numerical simulations of the plasma flow in the open field line tube are neces- 
sary to prove this idea. Given nonstationary generation of the secondary plasma 
the particles are confined to the separate clouds, and the fastest particles of a 
cloud can outstrip the slower particles of the previous cloud giving rise to the 
two-stream instability (Usov 1987, Ursov & Usov 1988). Up to now it is not 
clear whether the instability initiated in such a way can account for pulsar radio 
emission, since the nonstationarity of plasma generation is not studied in detail 
yet. The present paper suggests one more possibility of two-stream instability 
in pulsars, which is based on the selective energy loss of particles as a result of 
resonant inverse Compton scattering. Note that in this model the instability 
arises naturally, with no additional assumptions being involved. 

4 Gamma-ray luminosity provided by resonantly 
upscattered Compton photons 

The photons produced by the resonant inverse Compton scattering have the 
energies 

E - je B mc 2 ~ B 12 j2 McV, (14) 

so that typically E ~ 1 — 10 MeV. Note that the curvature gamma-photons 
produced in the polar gap as well as the photons upscattered by the primary 
particles have essentially higher energies, E ~ 100 MeV. Hence, the resonant 
scattering by secondary plasma results in an additional low-energy component 
in pulsar gamma-ray spectrum. The spectrum of this component in the case 
of some specific distribution functions of the scattering particles is obtained by 
Daughcrty & Harding (1989). In general, the low-energy tail of this component 
spreads even to the X-ray band on account of the photons resonantly scattered 
at distances z > R, where the magnetic field strength decreases essentially. 
However, these photons are very few, since at z > R the scattering rate decreases 
significantly due to the shrinking of the solid angle subtended by the neutron 
star. 

Given that the scattering is efficient, most of the energy of the secondary 
plasma should be transferred to the low-energy gamma-rays. The luminosity 
provided by the upscattered photons can be estimated as follows: 

L = nSmc 3 Aj, (15) 

where S is the cross-sectional area of the open field line tube, 

S=*-f, (10) 
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Rl is the light cylinder radius, A7 the difference between the Lorentz-factors 
of the particles, which mainly contribute to the final and initial plasma energy, 
A 7 ~ 10 2 - 10 3 . Substituting Eqs. (12) and (16) into Eq. (15) we find: 

L = 0.942 •10 3 ° gl2 ffi 373 CT gs/s. (17) 

Although in the particle rest frame the photons are scattered in all directions, 
in the laboratory frame they are beamed along the particle velocity. So the 
opening angle of the gamma-ray beam is given by ip = 3\/R/Rl- The averaged 
observed photon flux, F, is related to the luminosity as 

F " 2^<PE^E' (18) 

where VL — 7r<p 2 /4 is the solid angle occupied by the beam of upscattered pho- 
tons, d the distance to the pulsar, AE the energy band. Taking into account 
Eq. (17), Eq. (18) can be rewritten as 

5 /2 

F = 3 ■ lO-io gpf^|p photons/(cm 2 • s ■ kcV), (19) 

where d 3 = d/10 3 pc, E 6 = E/l McV. 

For most pulsars the flux given by Eq. (19) is too low to be detected. 
At present the detectors of the Compton gamma-ray observatory are the most 
sensitive to the low-energy gamma-ray emission (OSSE at 0.05-10 MeV and 
COMPTEL at 1-30 MeV). At 1 MeV the source sensitivity of OSSE is only 
2 • 10~ 7 photons/(cm 2 • s • keV) (Gehrels & Shrader 1996). In the observations 
reported by Kuiper et al. (1996) the flux from Geminga in the band of 3-10 
MeV was found to be 10 _7 E^~ 2 photons/(cm 2 • s • keV). Substituting Geminga 
parameters (P = 0.237 s, B i2 = 3.3, d 3 = 0.15) into Eq. (19) one can obtain the 
flux provided by the upscattered Compton photons: F = 3.8 • 10~ 7 n^j^R^ 2 E§ 2 
photons/(cm 2 • s • keV), which is consistent with the one detected. 



5 Conclusions 

We have investigated resonant inverse Compton scattering by secondary pulsar 
plasma. The process is found to cause the efficient energy loss of the secondary 
particles given the neutron star temperatures > 10 5 K, so that our results are 
applicable to most pulsars. For the resonant character of the scattering, the 
energy loss depends strongly on the initial particle energy. At 70 ~ 10 2 — 10 3 
the scattering is the most essential, the final Lorentz-factors of the particles 
being independent of the initial ones. 

The distribution function of the secondary plasma is significantly altered 
by the resonant inverse Compton scattering. It is shown that ultimately the 



9 



distribution function becomes two-humped. The main peak at 7 ~ 10 2 is very 
sharp. It is formed by particles which suffered severe energy loss on account 
of the scattering. Another hump is sufficiently broad. It is associated with the 
particles, whose Lorentz-factors are almost unaltered by the scattering. The 
two-humped distribution function of the plasma particles is known to be unsta- 
ble. It is shown that at pulsar conditions the two-stream instability develops 
readily and leads to an essential increase of plasma oscillations, which are likely 
to be transformed into radio emission. 

We have also estimated the gamma-ray flux provided by the upscattered 
Compton photons. The resonantly scattered photons appear to gain the energies 
of 1-10 MeV forming an additional low-energy component in pulsar gamma-ray 
spectrum. 
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Figure 1: Evolution of particle Lorentz- factor with the distance for initial 
Lorentz-factors, 7o = 100 (a) and 7o = 3000 (b); B 12 = 1, z Q /R = O.Ol.The 
solid lines are plotted for the case of scattering the photons from the whole 
neutron star surface with the temperature 5 • 10 5 K, whereas the dashed lines 
refer to the case when the contribution of the hot spot photons is also taken into 
account(here the hot spot radius is 5 • 10~ 3 i? and the temperature is 3 • 10 6 K). 
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Figure 5: Evolution of the waterbag distribution function on account of reso- 
nant inverse Compton scattering: a for various temperatures, T 6 ; Bi 2 = 1, b for 
various magnetic field strengths, Bi%\ Tg = 0.5; the initial distribution function 
is shown by the dashed line 
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Figure 6: The same as in Fig. 5 for the initial distribution function (5). 
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